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Abstract 

Theoretical analysis is presented on quantum state evolution of polarization light waves at frequencies Uo 
and bJe in a periodically poled nonlinear crystal (PPNC). It is shown that the variances of all the four Stokes 
parameters can be squeezed. 
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1 Introduction 

During the last decade much attention has been paid to the realization of quantum information protocols [1] 
such as quantum teleportation, quantum cryptography. These protocols are based on the methods of nonlinear 
quantum optics. Nonlinear optical sources [2] play an important role in the generation of nonclassical states 
of light [3] and realization of optical quantum information protocols. The nonlinear optical processes such as 
degenerate parametric down conversion (type I and II processes) [2,3] and the Kerr effect [2,3] are used to 
create nonclassical states of light (squeezed states, polarization squeezed states and entangled states). The 
variance of one of the four Stokes parameters of polarization squeezed light is less than the corresponding 
value for the coherent state. Traditionally the degenerate parametric process (type II) and the Kerr effect 
are responsible for the generation of polarization squeezed states or polarization entangled states in ordinary 
nonlinear crystals. One can achieve suppression of variances of at least one of the Stokes parameters [5-9] 
(5*0, Si, Sa) in the type II process by using ordinary nonlinear crystal. Most of the quantum information 
protocols are based on type II process [1] and Kerr effect, which are used to generate entangled states. The 
entangled states are used in the realization of quantum teleporation and quantum cryptography protocols. 
Experiments on quantum teleportation and quantum cryptography are performed by using ordinary nonlinear 
optical crystals with second and third order nonlinear susceptibilities. In the past few years, some experiments 
in the creation and realization of nonclassical and entangled states are performed by using PPNCs with second 
order nonlinear susceptibilities. The PPNCs ]10,11], which have many interesting advantages as compared to 
ordinary nonlinear crystals were proposed by Bloembergen and co-authors in 1962. The main advantages of 
PPNCs against ordinary nonlinear crystals are: the quasi-phase-matching condition between the interacting 
waves; the highest nonlinear susceptibility coefficient can be used; multi-mode interaction of optical waves. 

Recent experiments on quantum noise reduction [12,13] and the generation of entangled states J14-16] promis- 
ing the applications of PPNCs in the realization of optical quantum information protocols. These experiments 
were based on type I [12-14] and II [16] processes. It should be noted that the parametric down conversion 
(type I) and frequency sum generation processes have been studied theoretically Jsee for instance, J4] and the 
references therein] and experimentally [12-16] very well. The experiment on generation of time-bin and energy- 
bin entangled states using the parametric down conversion process (type I, i.e. 2tJe = u>e+ lOe) in a PPNC was 
demonstrated by Gisin and co-workers [14]. In this experiment authors claimed the higher energy conservation 
from the fundamental beam 2u]e to modes with frequencies uje and uje- In a PPNC, the possibilities of type II 
process is much more complicated as compared to ordinary nonlinear crystals. The following are the some of 
possible type II nonlinear processes which can be realized in a PPNC J4]: (a) o^o + i^e — 2L0e (one can achieve sup- 
pression of maximum three variances of Stokes parameters under certain conditions. The same can be achieved 
in an ordinary nonlinear crystal), (b) Uo + uje — 2uJe and 2ue + ujo = 3cJe (one can achieve suppression of all 
the four variances of Stokes parameters under certain conditions. The same can not be achieved in an ordinary 
nonlinear crystal), (c) uoo -\- uJe — 2a;e, 2LJe ~\- Uo ~ 3a;e, SctJ^ + uoe = 4ujo and 2uje + 2uje = 4:UJo (one can achieve 



higher suppression of all the four variances of Stokes parameters under certain conditions as compared to case 
(b). The same can not be achieved in an ordinary nonlinear crystal). The type II process (a) in PPNC can be 
accompanied by considering other nonlinear processes ((b) or (c)) [4]. The nonlinear process (a) differs with 
the same nonlinear process in an ordinary bulk crystal by a higher energy conversion rate of fundamental mode 
{2iUe) into degenerate (orthogonal) modes {ujo and tJe) [16]. The recent experiment on generation of polarization 
entangled states (type II process (a)) [16] demonstrated a high energy conversion rate from fundamental mode 
{2uje) into degenerate (orthogonal) modes {uJo and uie). The nonlinear processes (b) or (c) can be realized in 
a single PPNC but not in a single ordinary nonlinear crystal. All these nonlinear processes (b) or (c) can be 
quasi-phase-matched at certain coherent lengths [4]. Here we will study the generation of polarization squeezed 
states based on type II process (c) in PPNC with second order nonlinear susceptibility. 

The main goal of this work is to show that PPNCs can suppress all the four variances of Stokes parameters 
below the standard quantum limit. The theoretical work that we present here is the first, to the best of our 
knowledge, to propose PPNCs for the generation of polarization squeezed light. 

The structure of the paper is as follows. Section 2 describes the optical nonlinear processes and their 
Heisenberg equations of motions. Section 3 studies the behaviour of mean photon numbers of degenerate 
polarization (orthogonal) modes at frequencies uJo and uie. Section 4 analyzes the variances of Stokes parameters. 
The final section summarizes the results obtained in sections 3 and 4. 



2 Equations of motions 

We consider the five-frequency interaction of co-propagating light waves in a PPNC (see Fig.l). The four 
interaction processes of light waves at frequencies uJo, i^e, StJe, StUe, and Aluo are ]4] 

5ki = /c2e — kio - fcie + miGi = Afci + miGi, (1) 

LUo + 2ljJe = Side, 

5k2 = /C3e — fcio - k2e + m2G2 = Afc2 + m2G2, (2) 

UJe + 3LJe = 4lJo, 

(5^3 = fc4e - fcie - fcse + msGs = Afca -f m^Gs, (3) 

2uJc + 2uJe ~ 4cJo, 

5k4 = fc4e — 2k2e + ITliGi = Afc4 + UliGi, (4) 

where Afcj=i_2,3,4 is a phase mismatch for an ordinary nonlinear bulk crystal; /cj„ is a wave number of interacting 
modes (n = o stands for ordinary and n = e for extraordinary); Gj — 2-kAJ^ is the modulus of the vector of a 
reciprocal lattice with period Aj; mj = ±1, ±3, . . ., is the quasi-phase-matched order. 

The nonlinear process (1) describes the splitting of a photon of frequency 2uje into two photons of orthogonal 
polarizations with degenerate frequencies lOo and oje. The second process (2) describes the sum frequency 
generation, i.e. photon of frequency 2uje combines with the photon of frequency ujo, which gives rise to a photon 
of frequency SoJe. The third (3) and fourth (4) processes are responsible for the generation of the fourth harmonic 
with two different ways, i.e. a photon of frequency uJe and a photon of frequency SuJe combine, and a photon 
with frequency Aujo appears as a result or the same can be achieved by the combination of two photons with 
frequencies 2uJe. 

It has been shown that the nonlinear processes (1-4) can be simultaneously quasi-phase-matched ]4] in a 
single domain structure (Gi = G2 = G3 = G4) or at certain coherent lengths i^^^. The processes (1-4) under 
study can be simultaneously quasi-phase-matched with the following relationship 

Llol = 9i'l. (5) 
The processes (1-4) can be described by the following interaction Hamiltonian 

Hi[z) = ng(z)[^iaioaiea2^e + ^2aioa2ea3^e + iia^eai^a^^e ^ + C4a2<i"4oe + HC\, (6) 

where h is Planck's constant, ajn{a^^) is the annihilation (creation) operator of photon of jnth mode at frequency 
jiOn] £,j is the nonlinear coupling coefficient; g{z) is the periodic function equal to +1 or —1 at the domain 
thickness I — A/2; HC denotes Hermitian conjugate. The operators ajn{a'j'^) obey the following commutation 
rules 



[aj„,ap„] = Sjn^pn, with j,p= 1,2,3,4 



(7) 



The interaction Hamiltonian (6) can be averaged over the period A, if the interaction length z is much more 
than the period of modulation A, i.e. tv. Then the interaction Hamiltonian (6) takes the form 

Hi = h[-yiaioaiea2e + 72aioa2ea3e + TsOieOaeoJo + -yialeato + HC], (8) 

where 



+A/2 

7j = ^ J g{z)exp{±iAkjz) = 2^j/(Timj) 



-A/2 



The Heisenberg operator equations corresponding to the interaction Hamiltonian (8) are given by 

= 72a2el3e +7iaiea2e, 

dz 

.da\e , + , , ^ + ~ 

dz 

.da2e „ . . , , 

t—. = 274a4oa2e + 72<IloCl3e + Jiaieaio, 

dz 

.doL'ie - . , . + . 

— = 72aioa2e + 73aiea4o, 

,rf04o -2 , - - /„N 

= 74a2e + 73aieCl3c- (9) 

Assuming the pump modes at frequencies uj2e., ^jJao are classical and non-depleted at the input of a PPNC, i.e. 

a.2e = A2e, (10) 
CLAo = Aio, (11) 

where A2eAo = e"^''^ are the complex amplitudes of the pump modes, we obtain the following linear system of 
equations, after the substitution of quantities (10) and (11) into (9): 

daio . , 

= -72a3e +7iaie. 

=73a3e +7iaio> 

= 72010 + 73aie- (12) 
For simplicity, we introduce a normalized interaction length parameter Q 

C = 27i- (13) 

The quantity (13) is introduced into the set of equations (12), which reduce after straightforward algebra to the 
set of equations 

daio , „ , 

— — = -fc2a3e + Hie, 

dC, 

= feflio + fcia^'g. (14) 

where ki = 73/71 and k2 = 72/71. The set of linear equations (14) is solved by applying the Laplace transfor- 
mation: 

aio(C) = Aii(C)aio -I- Ai2(C)aie + ^13(^)036, 

aie(C) = A2l(C)ai + A22(C)aie + A23(C)a3"e, 

a3e(C) = A3l(C)Olo -I- A32(C)aie + A33(C)a3e. (15) 



where (0) and 

q = y/l + kj - kl 

Aii(C) = i(— fci coshgC + cosh + ^i), 

Ai2(C) = ^(-fcife coshgC + fcife + gsinhqC), 

AisfC) = -^(— ijfe sinhgi^ + fci coshg^ — ^^^i), 

A21 (C) = —^{1 sinh — fci + fci A:2 cosh qQ , 
<1 

A22 (C) = -T (^2 cosh qC + q^ cosh qC, ~ kl), 

q 

A23 (C) = '^il^^ sinh (jC — cosh qC + ) 1 
A31 ((^) = -^{ik'z sinh — fci + fci cosh g^) , 

A32(C) = l(.fc.sinh,C + fc.cosh,C-M, 
1 2 

A33(C) = ^(l + g cosh - cosh gC). (16) 

If fcl = fc2 = 0, then the solution (15) corresponds to the conventional degenerate parametric down conversion 
process (type II) [3]. Using (15) we can calculate the statistical properties of interacting modes with frequencies 
ujo, uJe, 3uje Sit normalized interaction length 



3 Evolution of mean photon numbers 

The evolution of mean photon number of degenerate polarization modes at frequencies u>o and cJe in a PPNC 
are calculated by using (15) for the following initial conditions at the input of the PPNC: the two orthogonal 
modes are in polarized coherent states \aio >, |aie >, Q!io,e = |Q!io,e|e"^^°'= , where |Qio,ie| = \/ < A''io,e(0) > 
and the wave at frequency SuJe is in the vacuum state |0 > 

< Nio.eiO > = < aio\ < aiel < 0|a^o^^(C)aio,e(C)|0 > jaie > |aio > (17) 

The expressions for mean photons (17) are 

< iVlo(C) >= A?2(C) + A?2(C)laie|' + 2All(C)Al2(C) COS {,Plo + 4>le) + \\AC)Wio\\ (18) 

< iVi.(C) AL(C) + A23(C) + AL(C)l«ie|' + 2A2i(C)A22(C) COS (0io + 01,) + A^i(C)|Qio|'. (19) 

It is well known that the parametric down conversion process (type II) depends upon the initial phases 0io,e of 
the polarized (orthogonal) coherent states. The values of mean photon numbers (18) and (19) depend upon the 
sum of initial phases, i.e. 4>io + 4>ie.- Under the condition 4>io + "i^ie = 2-ks, the mean photon numbers (18) and 
(19) start increasing rapidly with the growth of interaction length C,, under the condition <^io + 4>ie = tts, they 
start decreasing and then monotnically increasing (see Fig. 2) and s — ±1, ±2, .... 

Fig. 2 demonstrates the dependence of mean photon numbers < N-io > and < A^ie > on the normalized 
interaction length under the condition (^lo + 0ie ~ tt. It is seen that the behaviour of mean photon numbers at 
frequencies lOo (curve 1) and cJe (curve 2) are quiet different from the evolution of the mean photon number for 
the case of parametric down conversion (type II) in an ordinary nonlinear crystal. The mean photon number 
(curve 3) for the case of an ordinary nonlinear crystal is calculated by putting ki — kz — Q into the expressions 
(18) and (19) under the same initial conditions. 

The difference between the evolutions of mean photon numbers (18) (curve 1) and (19) (curve 2) illustrated 
in Fig. 2 is related with the interaction process (2), which is responsible for the sum harmonic generation, 
i.e, photon of frequency lOo combines with the photon of pump frequency 2uje, which gives rise to photon with 
frequency 3a;e. This can be easily seen by evaluating the mean photon numbers < NadC) > using (15) under 
the same initial conditions 

< iV3e(C) ^UO + A^2(C)l«le|' + 2A3l(C)A32(C) COS {c^lo + M + A^i (C) | ^lo | ' . (20) 

The curve (4) of Fig. 2 illustrates the growth of the mean photon number < N-ie > as the normalized interaction 
length increases. 



At earlier stages (normalized interaction length ^ ~ -r- 0.4), the mean photon number of frequency tJe 
decreases (see Fig. 2, (curve 2)) and later starts increasing monotonically. The decreasing of mean photon 
number < A''ie > is connected with the nonlinear process (3), which is responsible for the sum harmonic 
generation at frequency 4ljo and is realized by the combination of photons with frequencies uje and ScJe. The 
fourth harmonic generation process (3) is complicated as compared to (1) and (2). So, the interaction process 
(3) starts acting at later stages of interaction as compared to the interaction process (2). 

The calculations of photon variances < ANi^ ,,{Q > of polarization modes with frequencies lOo and iOe have 
shown that they are super-Poissonian [3]. The same can be seen in the ordinary parametric down conversion 
process (type II) in ordinary nonlinear crystals by putting the nonlinear coupling coefficients fci, ^2 equal to 
in the expressions (16) and substituting a3(0) = into the solution (15). 



4 Stokes parameters 

The polarization properties of orthogonal modes at frequencies loo and iOe can be analyzed by the Stokes 
parameters [8] 

S2(0 = ai,(C)aie(C) + a^e(C)aio(0, 

SsiO = i[a+ (C)aio(C) - &UO&ieiOl 



The Stokes parameters (21) obey the commutation relations of the SU(2) algebra [8] 

[5o(C),Si,2,3(C)] = 0; [5i(C),52(C)l =2iS3(C); 
[52(0, S3(C)] = 2iSi(C); [S'3(C),Si(C)] = 2iS2(C). 

The Heisenberg uncertainty relation for the Stokes parameters (21) is given by [8] 

< ASfiO X ASUO >> I < Sk > f,{i,j,k = l,2,3)(i / k). 



(21) 



(22) 



5o,i(C) = 

^2,3(0 = 

where — 



PO±(C) +Pl±(C)«3e«3e + P2±{C){"'lea3e + CL^Ate} + P3± (C)"i"le + P4± (C) {^100^ + ^lo"-:ie} 

+P5±(C){aioaie + a^^a^^} + P6±(C)a]*oaio, (23) 

'- i°'^[<lo{0{0-le ± alt} + gi(C){ai«3e ± aieffli} + I2{0{ate ± } + ?3(C){aioa3e ± CttoCltj 

+q4{0{aioate ± ato&ie} + g5(C){"L ± a?o }], (24) 

1, = i stand for S2{(), S'3(C) and 

Po±(C) = A?2(C)±Aii(C)±A^3(C), 
Pi±(C) = A?3(C)±Ai3(C), 

P2±(C) = Al2(C)Al3(C) ± A22(C)A23(C), 

P3±(C) = a?2(C)±aL(C), 

An(C)Ai3(C)±A2i(C)A23(C), 

All(C)Al2(C)±A2l(C)A22(C), 

P6±(C) = a?i(c)±aL(C), 

go(C) = Ai3(C)A23(C), 
Ai3(C)A22(C)+Ai2(C)A23(C), 
g2(C)=Ai2(0A22(C), 
Ai3(C)A2i(C)+Aii«)A23(C), 
Al2(C)A2l(C)+All(0A22(C), 

g5(C)=An(C)A2i(C). (25) 



P4±(C) 
P5±(C) 



ai(C) 

93(0 
94(0 



Further for 
expressions 
expressions 
(17)-(20) 



simplicity and clearness, we will write Pj (C)(j=o,i,2,3,4,5,6) as pj and 9^ (C)(j=o,i,2,3,4,5) as qj. The 
for the variances of the Stokes parameters (23) and (24) become very lengthy. So, we write down the 
of variances for the Stokes parameter Sj (0 under the same initial conditions applied for expressions 



< A5f (0 >=< S|(0 > - < S,(0 (j 



0,1,2,3) 



(26) 



The variances of Stokes parameters (26) are normalized by dividing them by the variance of the Stoke's parameter 
£•0(0), i.e. 



VAO - ^^m^, (27) 

where < AS'o(O) >= |aiop + For the case of coherent light, the relative variances (27) take values 1. If 

at least one of the relative variances (27) takes value less than 1, the light is said polarization squeezed. The 
relative variances V^(C) read 

Vb,l(C) = , 2,^, bo± +P2± +P5± + {P2± +2po±P3± +P3± +fl5±)|ai<i|^ 

< ASiiO) > 

+P3±\ctlef + 2(p2±P4± + 2po±P5± + P3±P5± + P5±P6±)\ctlo\\a.le\ COs{(l>io + (l>la) 
+ 4p3±P5±|Qlo||aie|^ cos (0lo + 0le) + 2p5± | Olo | ^ | Ole | ^ COs2((j!>lo + (f>i^) 
+ {P4± +P5± + 2po±P6± +P6±)l"lo|^ + 2(pi± + P3±P6±)\ceiof \ceief 
+4p5±P6±|aio|^|aie| cos (010 + 4>le) + pl±\aio\'^ 
-(PO±(C) + P3±{Q\aief + 2p5±(C)|"lo||aie| COs((/ilo + (/lie) + P6± (C) I "lo I ^ ) ^] • 

^2.3(0 = — TT^T— [2% + 2ql + ql + 2g| ± 2ql\ai,f cos (4<^ie) + (g? + 492^ + ql)\m,\^ + 2ql\ai,\* 

+ 2(gig3 + 25294 + 2g4g5)|aio||aiEl cos (0lo + 0le) + 4g2g4lQ!lo||Q!le|^ cos {(t>lo + ± 4g2g4lQ!lo||Q!le|^ COS (01O - 3(j!>le) 

+4q2g5!Qio|^|aiel^ COS {2(j)xo + 2(j)ie) ± 2(ql + 2q2q5)\aiof\aief cos (2<^io - 2(j!>ie) ± 45495 laio|^|aie| cos (3<?!)io - (f>ie) 
±2q'i\aio\'^ cos (4(j!>io) + 4g4<j5|cno|^|Qie| cos (0io + ^le) + (g| + gl + 4g5)|aio|^ + 2g||Q!io|^|aie|^ + 2q'i\aiof 
-{i°'^[q2\ai,f{e'^'^^' ± 6-'^'^'"} + g4|ai„||ai,|{e"^i°-"^i'= ± e-'^i°+'^i=} + g5|aiot'{e"'*^° ± e""'^!"}])"]. 

Under the conditions (pio + (pie = tt and ki < k2 expressions (28) and (29) become more favorable for the 
generation of polarization squeezed light and are evaluated for different initial mean photon numbers in polarized 
coherent states. The Fig. 3-5 demonstrate the evolution of expressions (28) and (29) as functions of the 
normalized interaction length Fig. 3 and Fig. 4 illustrate the relative variances of Stokes parameters 50,2 and 
5*1 and simulate the spontaneous parametric down conversion process, i.e. when the photons with frequencies 
LUo and iOe at the input of the PPNC are each having mean single photons |aio,eP = 1 in their orthogonal 
coherent states. Fig 3. shows that the variances of the Stokes parameters 50,2 start with the sub-Poissonian [3] 
statistics and later ( > 0.45 become super-Poissonian. The evolution of the variances of the Stokes parameters 
5*2 and S3 are almost the same. So, that is why we are not demonstrating the variance of S3. Moreover, the 
variances of the Stokes parameters So, 2, 3 do not differ from the same variances of the Stokes parameters for an 
ordinary nonlinear crystal. The latter one can be seen by putting the values of nonlinear coupling coefficients 
fei, fc2 equal to in the expressions (16) and substituting a3e(0) = into (15). 

The variance of the Stoke's parameter Si is more interesting due to its sub-Poissonian statistics, which is 
shown in Fig. 4. It should be noted that the variance of Stoke's parameter Si for parametric down conversion 
process (type II) in an ordinary nonlinear crystal is super-Poissonian. The latter can be checked by putting the 
values of fci = ^2 = in the expressions (16) and substituting a3e(0) = into (15). So, the PPNC can suppress 
variances of the Stokes parameters So, 1,2 under certain initial values of phases (0io + (pie = tt) of polarized 
modes and nonlinear coupling coefficients (fci < ^2). 

The expressions (28) and (29) are calculated under the same initial conditions but with different mean 
photons |aio,eP ~ W'^ in each polarized mode. The Fig. 5 shows the evolution of variances of all the four 
Stokes parameters So, 1,2, 3. It is seen that all the variances of Stokes parameters have sub-Poissonian statistics. 
After analyzing the Fig. 3-5, it is clear that the larger the mean photon numbers in degenerate polarized 
modes, the more the squeezing in the variances of all the Stokes parameter's. Moreover, the variances of Stokes 
parameters S2,3 are ~ 0, i.e. the generation of almost Fock states. 



5 Conclusion 

In this paper we have investigated the problem of the generation of polarization squeezed light in PPNC with 
second order nonlinear susceptibility. In the case of classical and non-depleted modes at frequencies 2(jJe and 
4uJo, the exact solution of the Heisenberg equations of motions is obtained. The solution is exact quantum 
solution showing new possibilities for the generation of polarized squeezed states of light. We have used this 
solution to calculate the mean photon numbers of interacting modes and variances of Stokes parameters. It is 
found that the evolution of mean photon numbers of interacting modes with frequencies lJo and cJe differ with 
the evolution of mean photon numbers of the same interacting modes in ordinary nonlinear crystals (degenerate 



parametric process (type II)). The reasons are found, which explain the peculiarities of degenerate parametric 
down conversion (type II) in PPNCs. 

Also, it is shown that all the four variances of the Stokes parameters can be squeezed (variances of the Stokes 
parameters are smaller than the value of variance of coherent state) simultaneously. Optimal initial conditions 
are found under which the squeezing of variances of the Stokes parameters can be obtained. 

So, PPNCs can be good candidates for the generation of polarization squeezed states of light, which can be 
used for the realization of quantum information protocols. 
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Fig. 1 Sektch of generation of polarization squeezed light in PPNC with second order nonlinear susceptibility. 
The light waves involved are described by their frequencies ujo, ^a, 2tt;e 
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Fig. 2 Mean photon numbers < iVio(C) > (1), < Nie{C) > (2), < N:ie{C) > (4) with frequencies ujo, t<Je, and 
SoJe, respectively and < Nio,e(() > (fci = = 0) (3), as functions of the normalized interaction length (. The 
curves (1-4) are calculated corresponding to the initial photon numbers < A'^io,e(0) 1 and < NseiO) >= 

3| 1 1 1 1 1 1 1 1 
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Fig. 3 Normalized variances Vb(C) (0) and V2(C) (2) of the Stokes parameters S'o(C) and ^(C)- Curves (0 
and 2) are calculated corresponding to the initial photon numbers < A'^io,e(0) >= 1 and < A'^3e(0) >= 
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Fig. 4 Normalized variance Vi{() of the Stoke's parameter Si{(). Curve is calculated corresponding to the 
initial photon number < A'^io,e(0) 1 and < A'^3e(0) >= 
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Fig. 5 Normalized variances Vo(C) (0) and Vi{Q (1), and V2(C) (2), of the Stokes parameters 5'o(C), •S'i(C) 
and 52 (C)- Curves (0-2) are calculated corresponding to the initial photon numbers < Nio,e{0) >= 10"^ and 
< N:i,{0) 



